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Abstract 

The mean spectrum and burst statistics of gravitational waves produced by a cosmological 
population of cosmic string loops are estimated using analytic approximations, calibrated with 
H-j ' earlier simulations. Formulas are derived showing the dependence of observables on the string 

£vq ■ tension Gu, in the regime where newly-formed loops are relatively large, not very much smaller 

, than the horizon. Large loops form earlier, are more abundant, and generate a more intense 

> : 

stochastic background and more frequent bursts than assumed in earlier background estimates, 
enabling experiments to probe lighter cosmic strings of interest to string theory. Predictions are 

in ; 

compared with instrument noise from current and future experiments, and with confusion noise 
from known astrophysical gravitational wave sources such as stellar and massive black hole binaries. 
Qh' In these large-loop models, current data from millisecond pulsar timing already suggests that Gfi 

is less than about 10 -10 , close to the minimum value where bursts might be detected by Advanced 

-4— > ■ 

c/3 ■ 

LIGO, and a typical value expected in strings from brane inflation. Because of confusion noise 

. ^ ' expected from massive black hole binaries, pulsar techniques will not be able to go below about 

V ■ 

$_i ■ Gji « 10 -11 . LISA will be sensitive to stochastic backgrounds created by strings as light as 



Gu ~ 10~ 15 , at frequencies where it is limited by confusion noise of Galactic stellar populations; 
however, for those lightest detectable strings, bursts are rarely detectable. For G\i > 10 , the 
stochastic background from strings dominates the LISA noise by a large factor, and burst events 
may also be detectable by LISA, allowing detailed study of loop behavior. Astrophysical confusion 
might be low enough at 0.1 to 1 Hz to eventually reach Gfi ~ 10 -20 with future interferometer 
technology. 

PACS numbers: 98.80.-k, 98.70.Vc,11.27.+d,04.30.Db 
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I. INTRODUCTION 



Cosmic strings have been studied for many years as a possible new form of mass-energy 
with new and distinctive astrophysical effects Originally conceived in field the- 
ory as defects caused by broken U(l) symmetries in Yang-Mills theories, they have recently 
re-emerged as possible structures of fundamental st ring theory, called cosmic superstrings, 
F-strings and D-strings Bfl S H, fl Q, Q, Q, H Q that tend to arise naturally in 
models of brane inflation [17[ . In both cases, a cosmic network of macroscopically extended, 
quasi-stable strings naturally forms via the Kibble mechanism, a process of quenching after 
cosmological inflation. Although strings are now known not to contribute significantly to 
the formation of cosmic structure, other gravitational interactions of even very light strings 
may be observable and give direct data on their fundamental nature. The most novel and 
conspicuous effect comes about from their principal decay process: the primordial network 
of long strings continually intersects with itself, spawning isolated, oscillating loops that 
ultimately radiate almost all of their energy into gravitational waves 

This paper estimates the spectrum and burst statistics of the cosmic gravitational wave 
background produced by loops of light cosmic string, and compares those predictions with ca- 
pabilities of current and planned experiments that span almost a factor of 10 12 in frequency: 
millisecond pulsar timing, the space interferometer LISA, and ground-based interferometers 
such as LIGO, VIRGO, and GEO. The main new results arise from renewed attention to the 
effects of the size spectrum of string loops: if a non-negligible fraction of loops stabilize when 
they are relatively large compared to the horizon, the surviving population forms earlier and 
loops are more abundant. As a result, the strings generate more gravitational waves, so that 
even very light strings are detectable. 

Although the fundamental physics differs widely for different types of strings, just two 
main parameters, depending on that physics, control their astrophysical behavior and grav- 
itational effects: the dimensionless mass per length or tension G[i (in Planck units with 
G = nipf anck , and setting the speed of light to unity) and the interchange probability p. The 
main aim here is to explore the sensitivity and relative strengths of the various experiments 
at small values of Gfi. The predictions here for both the noise sources and the string back- 
grounds are uncertain by at least a factor of a few, but are good enough to estimate the 
experimental sensitivity to (and currently allowed range for) string parameters, to identify 
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promising approaches for improvement, and to assess in general terms the science impact of 
the experimental projects. 

In addition, other parameters retained here, 7 and a, characterize uncertainty or con- 
troversy at present of certain aspects of the string behavior. The dimensionless decay rate 
7 is generally taken to be around 50 and is unlikely to differ from this by a large fac- 
tor. However, dynamical simulations of network behavior disagree about a, the size of 
newly formed loops in units of the Hubble scale H~\ and this strongly affects the predicted 
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311 ]. Recent work [27] suggests that 



that newly formed loops are smaller than the horizon by a factor of the order of a ~ 0.1, 
in contrast to earlier estimates of the order of a ~ G/i or even smaller — a difference of 
many orders of magnitude. The approximations adopted here are valid over a range of rel- 
atively large a, 0.1 > a > \0 5, ~yGfi, that is consistent with the new simulations designed to 
accurately estimate a (27], but that has not been included in recent background estimates. 
Within that range, the predicted spectra simply scale as flux oc a 1 / 2 . As shown below, larger 
a leads to more loops, stronger predicted backgrounds and hence greater sensitivity to light 
strings than earlier estimates. 

The analytical model of the network behavior, loop formation and loop radiation adopted 
here shares the basic framework of previous models. The most important difference is 
the choice of a large value of a, which leads to the approximation that loops dominating 
gravitational wave production form during the radiation era instead of the matter-dominated 
era, and because of the increased intensity, leads to consideration of lower Gfi than earlier 
studies. Some new features need to be added in turn because of the very low G/i, since 
some experiments see the fundamental frequencies (/ = 2/L, where L is the length) of 
loops that have a lifetime larger than the Hubble time today. Aside from those differences, 
the model here follows standard "one scale" approximations for string behavior. We adopt 



earlier calculations 
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as a normalization in high frequency limit, and scale these results 
to lower frequencies and lower G\x. 

The rate of observable beamed bursts of radiation from strin g cu sps is also estimated, 



32, 



33, 



34, 



351 ] . However, the 



and roughly agrees with extrapolations of other recent studies 
abundant loops arising from large a lead to a situation where bursts are also accompanied 
by a more intense stochastic background. The larger overall flux of stochastic radiation, 
especially from high redshift loops, is detectable at very low G/i, and for the lightest strings 
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detectable via their stochastic background, bursts are rarely detectable, both because the 
confusion is greater and because the fundamental frequencies are higher. (The detected 
frequencies of the stochastic background can be close to the fundamental frequencies of the 
most abundant loops, where they are not beamed substantially at all.) On the other hand 
somewhat heavier strings produce both an easily detectable stochastic background and a 
possibly detectable rate of bursts, detectable at high frequencies, especially with Advanced 
LIGO. In this case the unique identifiable signature of the bursts provides independent 
information about the loop parameters and behavior, and a validation of the source model. 

Another major difference from previous work on string detection is a more comprehensive 
comparison with sources of astrophysical confusion noise limits from populations of binaries. 
At low frequencies, including those accessible to LISA and to pulsar timing techniques, these 
will eventually be more important and fundamental than instrument noise in setting the 
ultimate attainable sensitivity. The confusion noise and instrument noise models for LISA 
frequencies are largely taken fromjsfil]. 

Current models of string-based inflation provide some context for the scale of string pa- 
rameters. Inflation is an effective theory, not a fundamental one. It describes the behavior 
and cosmological gravitational effects of a scalar inflaton field in terms of an effective poten- 
tial; the goal has always been to relate that field and potential, at a deeper level, to other 
fields including the Standard Model. In string theory, the recent worked examples with such 
a physical basis are models of brane inflation. An effective inflaton can be realized where the 
scalar field represents the separation of a pair of D3/D3 branes^J]. At the end of inflation, 
the branes annihilate and break U(l)xU(l) symmetries giving rise to cosmic Dirichlet strings 
(Dl-branes) and fundamental (F-) strings. In these models of brane inflation, if parameters 
are chosen to create the observed cosmic perturbations (5(f) ~ 10 -5 ) via the usual scalarquan- 
tum fluctuations, a wide range of striogteosioos is possible, HT* < G „ < HrilBQ, 
with a typical value Gfi ~ 10~ 10 . The range of p is not known but a variety of models 
yields a typical range 10~ 3 < p < The analysis here shows that current and planned 

gravitational-wave experiments probe this entire range of parameters and well beyond. In- 
deed, apart from the remaining doubts about a, the arguments above suggest that the typical 
value of Gfj, ~ 10~ 10 is already in conflict with pulsar timing data. From this perspective, 
current and future gravitational wave experiments contribute unique and useful information 
about inflation and its relation to new, stringy physics. 
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II. MODEL OF STRING LOOP POPULATIONS AND RADIATION SPECTRUM 



Many aspects of the behavior of strings in a cosmological context are well established 
The Kibble mechanism lays down local gauge strings initially as a set 
of nearly random walks with a very small comoving coherence length. These then evolve 
obeying the Nambu-Goto action; in an appropriate gauge, the local equation of motion of 
the strings is just the relativistic wave equation. On scales larger than the Hubble length 
H -1 , the network is "frozen" and stretches with the cosmic expansion. On smaller scales, the 
strings move with velocities of order c. They occasionally cross and exchange partners, with 
probability p at each crossing. In this way, closed loops separate from the network of long 
strings. It is this population of loops that dominates the production of gravitational waves. 
Each loop oscillates at a set of discrete set of frequencies 2n/L depending on its length L, 
and loses energy slowly due to gravitational radiation. We assume here that gravitational 
radiation is the main energy loss, which is the case for many varieties of fundamental strings. 

We adopt the familiar "one-scale" model of the loop population 20, 



21 



22 . The 



typical size of a newly born loop is parameterized conventionally as a/H. The fraction of to- 
tal horizon mass in such loops is of order G/i, times a numerical factor of order unity we take 
from numerical calculations of the scale-free, high redshift solution. The gravitational wave 
background at each frequency is a sum of contributions from loops at all redshifts contribut- 
ing at that frequency; here we approximate this by estimating the redshifted radiation from 
the brightest loops at each redshift. As we see below, the loop contributions fall naturally 
into two main populations: a high-redshift H population of long-decayed loops, radiating 
from close to the fundamental mode, and redshifted by a large factor; and a present-day P 
population, where the redshift is of order unity or less. The H loops that decayed at high 
redshift dominate the integrated flux at high frequency, in the form of a nearly Gaussian 
stochastic background. The P loops contribute significantly near their fundamental fre- 
quency, as well as a power law spectrum oc f~ 1 ^ above it; they dominate the low frequency 
background, and are also the source of rare brief burst events at high frequency. (Although 
this nomenclature is useful for discussion, one should bear in mind that the H and P con- 
tributions are comparable, within a factor of a few in total flux, over a range of about a 
thousand in frequency.) 

A population of loops of a variety of sizes and redshifts produces a background with a 
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nearly continuous spectrum. At any given time, the power in a broad band peaks near a 
frequency around the fundamental frequency of loops whose decay rate is H, 

f peak = 2H/ 1 Gfx (1) 



where 7 = 5O750 represents the dimensionless gravitational radiation power in all modes 2l( , 
estimated to be 7 50 ~ 1. At present, the peak of the spectrum is at about 

f P eak « 10- 10 (75oG/i/10- 9 )- 1 Hz. (2) 

At frequencies near f pea k, the background comes mainly from the present-day population P. 

We assume here that loops formed at high redshift, during the radiation era. This is 
true at all the frequencies we are concerned with, as long as the P population loops formed 
before matter-radiation equality at t eq , which requires: 

'jGfx/Ho < a/H eq , (3) 

or inserting numbers from concordance cosmology (see below), 

a > 10 5 7G>- (4) 

Thus the approximations used here can only be applied consistently to light strings and 
relatively large a. 

The intensity of the background in this situation is determined by the production rate 
of loops in the radiation era. Indeed, at sufficiently high frequencies, the observed radiation 
comes from loops which both form and decay in the radiation era. During the period between 
the formation of a loop and the time it decays, its energy does not redshift; thus the ratio 
of gravitational wave energy to radiation energy density is of order Gfi(Zf or . m /Z decay )p~ 1 , 
where Zf orm and Zd eC ay denote the inverse scale factors for formation and decay of typical 
loops, and the mass of newly formed loops is about a fraction Gfi of the total horizon mass. 



1Q: 



The spectrum in this limit is flat in energy density, and does not depend on frequency 



V GH = bsoGWV^L = lO-^G^/lO-^^VSso^ai 2 , (5) 

where Qgh refers to the energy density in gravitational waves per log / in the high frequency 
limit, in units of the critical density, and VLr denotes the total energy density in relativistic 
species. The coefficient Fl is a numerical factor of order unity characterizing the creation 



rate of new loop mass per Hubble time in a scaling solution; its value only depends on the 
mean equation of state at loop formation being that of relativistic particles (w = 1/3). The 
factor used here is taken so that the spectrum normalization agrees with numerical estimates 

n 

of 21J in the high / limit, and this is the normalization for all the spectra estimated here. 
A family of these H spectra is shown in Fig. 1 for various values of Gfi. 

The scaling p~ x comes from a very simple argument: to achieve a scaling solution in 
a one-scale model, each long string must spawn loops a rate independent of p, therefore 
the number of long str ings it needs to collide with varies like p~ l . This scaling has been 
questioned in detail^, but we adopt it here as a simple placeholder. 

This flat spectrum holds approximately down to the frequency f eq of waves coming from 
loops that are decaying at the epoch of matter-radiation equality. In concordance cosmology 
we have Z eq = 1 + z eq pa Qm/Qr ~ 3300, and H t eq pa 10~ 5 ; thus 

fj fpeak = (HotegZeg)- 1 W 30. (6) 

For feq < f < fpeak, the spectrum from the loops decaying at Z decay (f) increases above 
Qgh- The typical loops decaying during matter domination, but before Z pa 2, during which 
Z oc t~ 2/3 , scale like Z form oc i~ 1/2 , Z decay oc t~ 2/3 , and / oc Z^ cay t^ cay . Thus they produce 
a spectrum of redshifted radiation tt G oc (Z f orm /Z decay ) oc (/ / feq)' 1 ^ 2 ■ With this scaling, 
the radiation from intermediate redshifts, though not negligible, is however subdominant 
compared to the high frequency tail from the P population loops at Z < 2, which goes like 
oc J -1 / 3 ; that is why it makes sense to conceptually separate the H and P populations for 
describing the effects. 

For this analysis we neglect features in the spectrum near f pea k due to cosmic acceleration 
since redshift of order 1, which are not neatly accommodated by simple scaling arguments. 
Neglecting the brief accelerating phase does not at all affect the estimate of the H population 
or feq<f< fztti contributions to the stochastic background, which are normalized to the 
microwave background flux. It does alter the spectrum near f pea k and reduce the contribution 
of the P population, since the older universe at given H means currently-decaying loops 
of size L decay are slightly larger (hence rarer, relative to the radiation density) than in the 
matter-dominated scaling. The density of nearby loops is pegged to the (model-independent) 

— 1 /2 

radiation density, times a factor Zf orm /Z decay] this factor is proportional to the L de ^ ay} since 
the accelerating universe is about twice as old as the matter-dominated model of the same 
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Hq, the frequency at L decay is about half as large, their mass density is smaller by about 
and their gravitational-wave energy flux is smaller by about 2~ 3//2 . Thus the numbers 
below overestimate the burst rate from nearby loops and the mean flux at / < f pea k from 



those loops, by a factor of a few. Similar factors contribute to some o: 
between previous estimates of burst rate, for example between [34] and 
approximation is closer to 



the differences 
the present 



3^] . The estimates of stochastic background flux above f pea k are 
little affected by this approximation. 

Since loops do not radiate below their fundamental frequency, at low frequencies / < 
fpeak the background comes from loops whose lifetime exceeds Hq 1 . The energy density 
spectrum is given by the mass density Vt^f^ 1 ) of loops of size L « 2/ _1 , which in turn is 
oc Zj orm oc L -1 / 2 oc J 1 / 2 , times the fraction of their rest mass they radiate in a Hubble time, 
which in turn is oc L^ 1 oc /: 

n G ^n L {f- l ){f/f peak )^fi\ (7) 

This regime was often not considered in earlier treatments because unless G\i is very small, 
these fundamental frequencies are too low to be observed directly. Although the normaliza- 
tion of this contribution is overestimated here due to neglect of recent acceleration, the slope 
should be accurate. The f2 G oc f 3 ^ 2 scaling applies only down to the frequency corresponding 
to loops formed at t eq , which have / ~ 2.5 x 10 _12 a(7 1 Hz. 
Summarizing the above contributions, the spectrum is: 

tt G (f) « n GH [l + V36(f/f P eak)- 1/3 }, f > m P eak (8) 

^GH[VS0(f/f P eak)' 1/2 + V30(f / fpeak)'^} , fpeak <f< ^fpeak 

n GH V30(f/f P eak) 3/2 , 2.5 x KT^Hz <f< f peak . 

Spectra in Fig. 1 show contributions from the various populations at each frequency. The 
radiation in these models is much more intense than estimates based on smaller values of a, 
where the present-day population is made of loops which formed during the matter era. As 
a result, the current and future sensitivity to G/i is correspondingly greater. 

We have adopted the conventional size parameter a for the one-scale model, but in reality 
not just the characteristic size, but also the detailed shape of the initial size spectrum 
matters. For such small Gfi, even if almost all loops shatter into tiny pieces quickly, a small 
fraction of loops that stabilizes close to the horizon scale can still come to dominate the loop 
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FIG. 1: Predicted gravitational wave backgrounds from strings, and noise sources. Broad band 
energy density is shown in units of the critical density for ho = 1, as a function of frequency, for 
a = 0.1. Noise levels are shown for current millisecond pulsar data (MSP), and projected LISA 
sensitivity in maximum resolution and Sagnac modes. Confusion noise is shown for massive black 
hole binaries (MBHB), the summed Galactic binary population including binary white dwarfs 
(UB+WUMa+GCWDB+CV), and extragalactic populations of white dwarfs (XGCWDB) and 
neutron stars (XGNSB). Radiation from loop populations at high redshift (H) and present-day (P) 
is shown, labled by the value of G[i. Dotted curves show the contributions of z > 1 loops where 
they are subdominant to the P contributions. Current (MSP) sensitivity is at about G\i ~ 10~ 10 , 
and LISA will reach to around G\x ~ 10 -15 . 

population. The fraction of horizon-size loops that needs to stabilize in order for radiation- 
era loops to dominate near f peak is only « (Gn) l l 2 (VL M / \f30£l R ) ps 600(G / u) 1 / 2 (times some 
numerical factors), which is a small number for the light strings we are contemplating. Unless 
a is relatively large, accurate estimates of the radiation background will require estimates 
of loop spectra with very large dynamic range. 
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III. BURST DETECTION AND SELF-CONFUSION NOISE 



At frequencies / much higher than the fundamental mode, the radiation is increasingly 
more highly beamed and episodic. The most important characteristic of loop radiation 
far above the fundamental frequency arises from the catastrophes called cusps. |37( It can 
happen that in the course of a typical loop orbit, at some instants, at single points on the 
world-sheet, the velocity of the string is formally the speed of light. The vectors describing 
the derivatives of the string trajectories can be mapped onto curves on a unit sphere for left- 
and right-moving modes, and cusps occur at some point whenever these intersect. Since 
curves "typically" intersect on a sphere, it is thought that the behavior happens reasonably 
often, of order once per orbit for a typical loop. 

The behavior of cusps obeys simple, universal scaling deriving just from the geometry 
of the events. The metric strain obeys a universal, power-law waveform that includes a 
power-law tail of energy up to high /. Approaching time St from the moment of a typical 
catastrophe, in units given by the fundamental mode of the loop, the metric strain amplitude 
due to radiation beamed from a cusp varies like h oc St 1 / 3 , and is beamed within an angle 
6 oc St 1 / 3 . (That is, if observed at angle 9, the cusp behavior is smoothed out at St < 6 3 .) 
The time-integrated spectrum at / is dominated by the waveform at St ~ Thus the 
integrated energy is oc h 2 f 2 9 2 St oc f^ 1 / 3 . 

These events are called "bursts". (Even though the amount of energy radiated close to 
a cusp actually decreases as St 1 / 3 , the power within the beam, in a particular direction, 
does increase like St~ 2 / 3 until the smoothing sets in.) The distinctive signature of a single 
burst stands out from other astrophysical sources, and can be distinguished if it has enough 
power relative to other noise sources, including the confusion noise of other loops and distant 
bursts. As discussed below, even with the optimistic assumptions made here, bursts are not 
in general the most easily detectable feature of very light strings; on the other hand, they 
can be a powerful diagnostic even if they are fairly rare. 

The number of loops contributing to the background from population H is so large that 
its contribution is for all practical purposes gaussian. The number of loops contributing from 
population P is also large near f vea k but decreases at higher frequency since the radiation 
from each loop is concentrated in angle and in time. Among the P loops, some are much 
closer than the mean so they can occasionally be bright enough to poke above the background 
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and produce a detectable burst. We want to estimate the rate of the detectable bursts. 



The rate of 
Vilenkin 



burst events was estimated in the important papers by Damour and 



32, I33I 3J]. A more detailed analysis by Siemens et al. |3j| yields lower rates 
due to different (and more realistic) cosmology and more demanding criteria for detection. 
The rates estimated here are not as reliable as Siemens et al. and err on the side of large 
burst rates. Even so, we find that bursts are not an important observable for the very 
lightest strings, because they are so light that they are observed rather close to their fun- 
damental mode. We make several simplifying assumptions, such that all the radiation in 
the high-frequency tail is from cusps and that there is one cusp per orbit. This makes no 
difference in estimating the strength of stochastic backgrounds, but it gives an optimistic 
estimate of the frequency of detectable bursts. Event rates simply scale directly with these 
numbers. 

The total number N of loops in the P population in a volume V is about 

^ PcTit^G{fpeak)V /gx 

M L (f peak ) K) 

where Mi(/ pea fc) « /jl/ f pea k is the typical mass of one of the loops. Adopting a volume 
V = Ati/3Hq and the spectrum as estimated above for Qcifpeak) yields an estimate of the 
total number of individual sources contributing near the peak of the spectrum, 

N(f peak ) « lO^G^/lO-^-^VSso 172 ^! 2 - (10) 

Note that in the same volume there are about 10 9 bright galaxies; for strings lighter than 
Gfi ~ 10~ 9 , loops are more numerous than galaxies. 

It is important for the following that the loops are distributed uniformly in space; we 
would get a very different answer for burst rates if they clustered in the halo of the Milky 
Way like dark matter. But since the center of mass of most of the loops moves rapidly 
through space, due to the "rocket effect" of beamed gravitational radiation, their large pe- 



culiar velocity likely prevents them from easily binding into galaxies 18]. Maximal beaming 
predicts that a loop has velocity of order unity by the time it decays (in spite of cosmological 
redshifting over a loop's lifetime, since most of the recoil momentum is imparted coherently 
over a decay time). As long as a typical loop has v > 10~ 2 ' 5 , more than the typical bind- 
ing velocity of galaxy systems, a seemingly safe but still unproven assumption, the bulk of 
them remain unclustered. (As an aside, note that if a loop does by chance become bound 
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its acceleration is not typically enough to allow it to leave a galaxy unless by chance the 
beaming is carefully controlled directionally to match the orbit. Thus a small population 
of loops does become bound to galaxies and clusters. For very light strings, there might be 
enough bound to our galaxy to affect the estimate of burst rates.) 

Near f pea k, the radiation from each loop is radiated in a broad beam pattern and more 
or less continuously. At higher frequencies, the power in the f" 1 ^ 3 tail is beamed in angle 
and concentrated in time around the cusps. The fraction of loops beaming in our direction 
is oc 9 2 oc /~ 2 / 3 , while the fraction of time radiating in a band around / is oc therefore 
the number of loops contributing to the background at any given time from population P 
at frequency / > f peak is about 

N(f)^N(f peak )(f/f peak )~ 5 / 3 . (11) 

In order to detect a burst with a matched filter with high confidence, the power of the 
burst must substantially exceed the noise power in the rest of the stochastic background 
during the burst. Since the burst only lasts at each frequency for a time t ~ the power 
exceeds the background by a factor of about SNR 2 where SNR is the signal-to-noise ratio 
in the metric strain; because they are rare events, we will require SNR > 5 for a detection, 
and more for a waveform study. For N(f) » 1, this happens only when a burst occurs 
unusually close to us, so that it stands out from the background. We estimate the mean 
time between burst events that stand out from this background. 

In the first instance, assume that the only background is from the same loop population 
that is producing the bursts, that is, made mostly of bursts from the P loops at distance 
~ Hq 1 . In order to be detectable above this background, a burst has to be within a distance 
D, where DH SNR~ 1 N(fy 1 ^ 2 . The total number of loops (including those not beamed 
at us) within this distance is then 

N D » SNR- 3 N(f peak )N(f)- 3 / 2 . (12) 

Each oscillation (that is, in a time the probability of a given loop sending a burst 

at us at frequency / at some point is about (/ '/ 1 f pe ak)~ 2 ^'' ! - The rate of observable events is 
then 

r-^Nofif/f^)- 5 / 3 , (13) 
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leading to a mean time between observable events, 

r P>H » SNR^NU^Y'^f/U^y^r 1 (14) 

or 

r P>H « 6 x 10 10 sec SWi^/lO"")" 19 /^^^^^ (15) 

where we require SNR 5 = SNR/5 to be at least of order unity for a significant detection, 
and significantly higher for a confirmation of the cuspy waveform shape. We can ask for 
example whether a burst might be seen by LISA from very light strings. Since the optimal 
frequency in this case / = 10~ 2 ' 5 Hz, we conclude that for very light strings near the detection 
limit, bursts are likely to be rare, at best occasionally detectable in an experiment lasting 
several years. 

Now consider the regime where the bursts have to compete with a dominant H back- 
ground. In that case, the detection distance D is smaller; it is multiplied by the square root 
of the P to H flux ratio, ~ [30 1 / 2 (///p e afc) _1//3 ] 1/ ' 2 - The time between detectable events is 
longer by a factor oc D~ 3 , so 

r P<H » 8 x 10- 2 SNR 3 N(f peak )^ 2 (f I 'f peak )-V 3 /-\ (16) 

or 

r P<H » 6 x 10 8 sec SNi^G^/lO" 12 )- 13 / 1 ^-^(f/lO-^Ez)-^ 3 ^ 12 ^}. (17) 

Thus at around G/i ~ 1CT 11 and above, burst events are common enough to be studied at 
LISA frequencies above their own confusion background. 

Now consider the instrument-noise limited case, which will apply at high frequencies such 
as those observed with LIGO. Again, since we require the power of the burst to exceed the 
noise power in the time interval we can specify the instrument noise in equivalent 
power units, Qi, the intensity of an isotropic stochastic background that would match the 
instrument noise during the integration interval ~ Similarly to the previous case we 
write 

t~i>h>p = [^gh/^iT Z/2 t p< h (18) 

which becomes 

t i>h >p « 10 12 sec SA^(^ 2 /10- 3 ^ 

(19) 

13 



where we have normalized to typical / and (very roug hly) Qj for Advanced LIGO[38j. (We 



have scaled from the multi-year sensitivity level of Qi ~ 0.5 x 1( 
is proportional to the inverse square root of integration time 



8 since in Qj units the noise 



36(). This estimate suggests 
that even for large loops, bursts will be observable with Advanced LIGO only for Gu 10~ 10 
or larger, in rough agreement with the more detailed analysis of Siemens et al. 35]. As we 
have just seen, the prospects for detecting bursts are somewhat better with LISA. 



Within the generous errors we have allowed for 



Durst rates in this exploratory discussion, 



the rate estimated here agrees approximately with and [3^ if we evaluate at the value of 
a for which loops are forming at t eq and allow for different cosmological and filter-threshold 
assumptions. The general conclusion is that the stochastic background, and not bursts, is 
the the best way to find and constrain light strings, but there is a still a window where 
bursts might be found and provide a valuable signature of string activity. 

IV. DETECTION SENSITIVITY AND NOISE SOURCES 

The primary detection techniques under consideration for measuring the stochastic back- 



ground are timing of millisecond pulsars. [391. 14C . 
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which is sensitive to gravitational 



waves with frequencies of order the inverse observation time (with the greatest sensitivity 
at longest periods, of order tens of years), and the spaceborne Laser Interferometer Space 
Antenna (LISA), currently under development, which is sensitive to waves in a broad band 
from 10" 2 to 10" 4 Hz.£| 

Apart from the burst detection, ground-based interferometers are not competitive; even 



advanced LIGO will at best reach a level 



Q G h 2 ps 0.5 x 10 8 , and will not be the best 



tool to detect stochastic backgrounds from very light strings. 

Figure 1 shows an estimate of sensitivity from a current pulsar study; the data have 
not yet been converted into a formal limit, but this gives an estimate of the currently 
attainable sensitivity level using the Pulsar Timing Array (PTA). A current benchmark 
of MSP sensitivity is the maximal monochromatic fit by LommenQ to residuals in the 
cleanest pulsar data, PSR B1855+09, which yielded Q G h 2 < 2 x 10~ 9 at / = 1.9 x 10~ 9 Hz. 
This corresponds to the maximum density of a gravitational wave of that single frequency 
allowed by the data, assuming the wave is oriented optimally to produce the largest possible 
signal in the direction of PSR B1855+09. Although this is not yet presented as a rigorous 
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upper limit it is shown in Fig. 1 as a guide to the sensitivity level possible with current 
timing accuracy and pulsar behavior. 

The view adopted here is that the actual limit from the best current data will be within a 
factor of a few of this estimate. This view contrasts with that of Damour and Vilenkin[34], 
who drew attention to conspicuous noise in some of the plotted data stream from this pulsar, 
and adopted an interpretation of the complete dataset that would allow significantly higher 
backgrounds. In this situation however, the noisy data came from one telescope system 
(the Green Bank Telescope) and much more noise-free data from another (Arecibo); the 
interpretation adopted here is the natural one based on only data from the best understood 
(and lowest noise) system. 

The sensitivity curve is also shown for LISA in two ways: the upper curve is the sensi- 
tivity to a stochastic background at the maximum frequency resolution (of the order of the 
inverse observation interval); the lower curve is the sensitivity, after three years, including 
information from a broad band, assuming that stronger confusion backgrounds are resolved 



into sources, and using the symmetric Sagnac combination 



44, 45. 



47| as a calibration of non- 



gravitational wave backgrounds. The Sagnac technique can use information from a broad 
band to measure a broad band background, and thereby gains up to two orders of magnitude 
in sensitivity for LISA to broadband stochastic backgrounds, depending on frequency 3(| . 
Depending on details of the mission still under development, the actual sensitivity below 
10~ 4 Hz may be significantly worse than shown. 

In addition to the instrumental noise, the detection of backgrounds is limited by astro- 
physical confusion noise. This ultimately limits the sensitivity of PTA, and of LISA at low 
frequencies. 

Consider for illustration a binary system losing its binding energy predominantly to grav- 
itational radiation. It radiates, over a decay time, the binding energy per mass E »s 2GM/ R 
in a broad band at the orbital frequency, / ~ (M/i? 3 ) 1 / 2 . Thus a homogeneous population 
of binaries over time creates a background j^S] with oc E oc f 2 ^. At low frequencies the 
background is confused because of a very large number of sources; at higher frequencies, the 
sources are rare and more intense, so that an experiment can separate out the individual 
sources by resolving them in frequency. The curves in Fig. 1 terminate at high / at around 
the frequency where the confusion from these backgrounds will be resolved out after a three 
year mission. 
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Backgrounds from known populations of binaries shown in Fig. 1 in the LISA band 
are taken fromj^; for convenience, the Appendix recapitulates assumptions and caveats in 
those estimates. At lower frequencies, the confusion background from massive black hole 
binaries 49, 50, 51, 5^] will ultimately limit the sensitivity of any pulsar experiment; indeed 
this background may be detected in the near future using the improved Pulsar Timing 
Array. That population of massive binaries will be resolved at higher frequency by LISA, 
and individual binaries will be detected with very high signal-to-noise ratio. In the lower 
part of the LISA band, confused signals from various populations of Galactic binary stars 
dominate the noise. In the upper part of the band, the Galactic binaries can be resolved 
out, but at these higher frequencies an irreducible background remains from white dwarf and 
other compact object binaries in other galaxies. At still higher frequencies of the order of 0.1 
to 1 Hz, above the LISA band but accessible in principle to a similar mission with shorter 
™fl Q, the white dwarf binaries disappear to „, and even the extragaiaetie 
neutron star and black hole backgrounds might eventually be resolved out. It could be that 
this band would allow many orders of magnitude improvement in sensitivity (although 
the possibility remains of other astrophysical noise from faint numerous sources that would 
limit the measurement.) 

It is worth noting that, if the string background is detected, it will only be possible to 
separate the very numerous loop sources in the rare bursts, so the stochastic loop background 
will be the dominant source of noise for LISA at whatever frequencies it is detected. 



V. CURRENT AND FUTURE DETECTABILITY 



Unless otherwise noted, for simplicity the following comments are made for a — 0.1, p — 1, 
as plotted in Fig. 1. At fixed /, the predicted loop radiation signal scales in this model in 
the three frequency regimes like: /i 1//2 a 1//2 p -1 (high /, H population dominated); /i 1 / 6 *^ 1 / 2 ^ 1 
(intermediate /, P population dominated); /x 2 a 1 / 2 p~ 1 (low /, P population dominated). 
These can be used to generalize to other combinations of /x, a, p. 

The current data on PSR B1855+09, with tt G h 2 < 2 x 10~ 9 at / = 1.9 x 10~ 9 Hz, appear 
to set an upper limit on Gfi of about 10~ 10 : at this value, the predicted flux lies about 
an order of magnitude above the monochromatic flux limit |43|| . At Gfi = 10~ 10 , the pulsar 
bound comes from close to (or just above) f pea k, so the background is truly stochastic, with 
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no detectable burst events expected. Note that at Gfi = 1CT 10 and above, the predicted flux 
at this frequency scales like f2c oc {Gii) l ^a l / 2 p~ l , so the limit on Gfi is quite sensitive to 
factors of order unity in the precise value of the limit on Qq. The limit also gets weaker for 
smaller a. If it is true, as assumed here, that a ~ 0.1, this limit is already in the middle of 
the range expected in brane inflation models, so a more careful evaluation of the bound is 
well motivated. That analysis should include a reasonable prior expectation of the predicted 
spectrum, which in this regime is not flat (and includes, in particular, features near f pea k 
due to recent cosmic acceleration.) 

At this level, bursts from the P loops may be detectable with Advanced LIGO data[35]. 
Such detections would serve both to confirm the stringy nature of the background source, 
and to break degeneracies in the parameters by separating H and P populations. 

Although it appears that a value of Gfi significantly above 10~ 10 is already ruled out, a 
different behavior sets in just below Gfi = 10~ 10 : since f pea k for Gfi < 10~ 10 ' 3 lies above 
the frequency where the pulsars are most sensitive, the turndown in the loop spectrum at 
/ < fpeak makes it significantly more difficult to set limits below Gfi ~ 10~ 10,3 . Thus, the 
current limit probably is within a factor of a few of Gfi = 10~ 10 . Future projects will achieve 
better sensitivity by improving timing residuals and increasing the number of pulsars. They 
will eventually be limited however by the MBHB background, which will limit sensitivity 
for significant detection or upper limits to not better than about Gfi ~ 10~ n . 

The smallest Gfi reachable by LISA appears to be in the Sagnac mode in a window 
around 10~ 4 Hz. At Gfi ~ 10~ 15 the background lies about a factor of three above the 
binary confusion flux, which in turn is about twice the broadband (Sagnac) instrument 
limit. Again, this is close to f pea k so it is stochastic and dominated by the P population, 
with little opportunity to see bursts. 

For higher Gfi than this, the loop background is a significant additional noise source for 
LISA. Although the P radiation dominates the total power in part of the LISA band for 
10~ 15 < Gfi < 10 -13 , as we have seen the individual burst events are seldom detectable as 
they are hidden in the confusion noise. For larger Gfi than about 10 -12 , the confusion noise 
across the entire main LISA band — about a factor of 100 in frequency — is dominated by 
stochastic H population radiation that is the first and most conspicuously detectable effect 
of the strings. For Gfi greater than about 10~ n , LISA-detectable bursts above this noise 
can occur often enough that their statistics can be used as a diagnostic of string properties. 



17 



34 



hese results confirm the generally robust prediction of burst events from strings [32L 133 . 
35| but also show that for small Gfi, the mean stochastic background is a much more 
sensitive probe than bursts. 

It is possible that confusion backgrounds are small enough to allow detection for Gfi < 
1CT 20 if a suitable detector could be built at 0.1 to 1 Hz , with sensitivity enough to resolve 
all the extragalactic binary sources and reach an instrument noise limit on the order of 
Q G h 2 ~ 10 -15 , the level required to search for inflationary graviton perturbations at zero 
spectral tilt (see e.g., 0, is of course also possible that unresolvable backgrounds 

not yet modeled could make this level unreachable.) By the same token, strings would have 
to be at least as light as this to avoid introducing a contaminating foreground that would 
make the inflationary waves undetectable. 

Acknowledgments 

This work was supported by NSF grant AST-0098557 at the University of Washington. 



T. W. B. Kibble, J. Phys. A 9, 1387 (1976) 

Y. B. Zeldovich, Mon. Not. Roy. Astron. Soc. 192, 663 (1980). 

A. Vilenkin, Phys. Rev. Lett. 46, 1169 (1981). 

A. Vilenkin, Phys. Lett. 107B, 47 (1981) 



M. B. Hindmarsh and T. W. B. Kibble, Rept. Prog. Phys. 58, 477 (1995), |hep-ph /9411342 
A. Vilenkin, E. P. S. Shellard, Cosmic Strings and Other Topological Defects, Cambridge 
University Press (2000) 

N. Jones, H. Stoica, and S. H. H. Tye, JHEP 07, 051 (2002), hep-th/0203163 

S. Sarangi and S. H. H. Tye, Phys. Lett. B536, 185 (2002), |hel>th/0204074[ 

N. T. Jones, H. Stoica, and S. H. H. Tye, Phys. Lett. B563, 6 (2003), |hep44i/0303269 

E. J. Copeland, R. C. Myers and J. Polchinski, JHEP 0406, 013 (2004) 

|arXiv:hep-th/0312067| . 

M. G. Jackson, N. T. Jones, and J. Polchinski, (2004), |hep-th /0405229 
J. Polchinski, 2004 Cargese Lectures, |arXiv:hep-t h/0412244 (2004) 



18 



J. Polchinski, arXiv:hep-th/0410082 



A. C. Davis, T. W. B. Kibble, Contemp. Phys. 46, 313 (2005) 
T. W. B. Kibble, |arXiv:astro-ph/0410073| 

A. Vilenkin, |arXiv:hep-th/0508135| 

S. Kachru, R. Kallosh, A. Linde, J. Maldacena, L. McAllister and S. P. Trivedi, JCAP 0310, 
013 (2003) |arXiv:hep-th/0308055| . 

C. J. Hogan and M. J. Rees, Nature 311, 109 (1984). 

T. Vachaspati and A. Vilenkin, Phys. Rev. D35, 1131 (1987). 

D. P. Bennett and F. R. Bouchet, Phys. Rev. D 41, 2408 (1990); 
R. R. Caldwell and B. Allen, Phys. Rev. D 45, 3447 (1992). 

R. R. Caldwell, R. A. Battye and E. P. S. Shellard, Phys. Rev. D 54, 7146 (1996) 
|arXiv:astro-ph/9607130"] . 

B. Allen and E. P. S. Shellard, Phys. Rev. Lett. 64, 119 (1990). 

D. Austin, E. J. Copeland, and T. W. B. Kibble, Phys. Rev. D48, 5594 (1993), 



hep-ph/9307325. 



X. Siemens, K. D. Olum and A. Vilenkin, Phys. Rev. D 66, 043501 (2002) 
|arXiv:gr-qc/0203006| . 

X. Siemens and K. D. Olum, Phys. Rev. D68, 085017 (2003), |g^qc/0307113 



V. Vanchurin, K. D. Olum and A. Vilenkin, arXiv:gr-qc/0511159 



C. Ringeval, M. Sakellariadou and F. Bouchet, arXiv:astro-ph/0511646 



C. J. A. Martins and E. P. S. Shellard, Phys. Rev. D 73, 043515 (2006) 
|arXiv:astro-ph/051l"792] . 

M. Sakellariadou, JCAP 0504, 003 (2005) |arXiv:hep-th/0410234| . 

A. Avgoustidis and E. P. S. Shellard, Phys. Rev. D 73, 041301 (2006) 
|arXiv:astro-ph/0512"582] . 

T. Damour and A. Vilenkin, Phys. Rev. Lett. 85, 3761 (2000), gr-qc/0004075 , 
T. Damour and A. Vilenkin, Phys. Rev. D64, 064008 (2001), |gr^qc/0104026| 
T. Damour and A. Vilenkin, "Gravitational radiation from cosmic (super)strings: Bursts, 
stochastic Phys. Rev. D 71, 063510 (2005) |arXiv:hep-th/0410222| . 
[35] X. Siemens, J. Creighton, I. Maor, S. R. Majumder, K. Cannon and J. Read, 
arXiv:gr-qc/0603115 



19 



C. J. Hogan and P. L. Bender, Phys. Rev. D 64, 062002 (2001) [arXiv:astro-ph/ ]. 

N. Turok, Nucl. Phys. B242, 520 (1984). 

C. Cutler and K. S. Thorne, |arXiv:gr-qc/ 0204090 

V. M. Kaspi, J. H. Taylor and M. F. Ryba, Astrophys. J. 428, 713 (1994). 
S. E. Thorsett and R. J. Dewey, Phys. Rev. D 53, 3468 (1996). 

M. P. McHugh, G. Zalamansky, F. Vernotte and E. Lantz, Phys. Rev. D 54, 5993 (1996). 
A. N. Lommen and D. C. Backer, Astrophys. J. 562, 297 (2001) arXiv:astro-ph/0107470| . 
A. N. Lommen, |astro -ph/0208572 , 

J. W. Armstrong, F. B. Estabrook, & M. Tinto, ApJ 527, 814 (1999) 
M. Tinto & J. Armstrong, Phys. Rev. D 59, 102003 (1999) 

F. B. Estabrook, M. Tinto & J. W. Armstrong, Phys. Rev. D 62, 042002 (2000) 

M. Tinto, J. W. Armstrong & F. B. Estabrook, Phys. Rev. D, 63, 021101 (R) (2000) 

E. S. Phinney, |a7Xiv:astro-ph7oT0 8028 . 

M. Rajagopal and R. W. Romani, Astrophys. J. 446, 543 (1995) arXiv:astn>ph7 9412038 . 



A. H. Jaffe and D. C. Backer, Astrophys. J. 583, 616 (2003) |arXiv:astro- ph/0210148 . 



J. S. B. Wyithe and A. Loeb, Astrophys. J. 590, 691 (2003) |arXiv:astro-p h/0211556 . 

M. Enoki, K. T. Inoue, M. Nagashima and N. Sugiyama, Astrophys. J. 615, 19 (2004) 

|arXiv:astro-ph/0404389"] . 

N. Seto, Phys. Rev. D 73, 063001 (2006) arXiv:gr-qc/0510067| . 

J. Crowder and N. J. Cornish, Phys. Rev. D 72, 083005 (2005) |arXiv:gr-qc/b 506015 . 

T. A. Prince, M. Tinto, . L. Larson and . W. Armstrong, Phys. Rev. D 66, 122002 (2002) 

|arXiv:gr-qc/0 209039 . 

A. Krolak, M. Tinto and M. Vallisneri, Phys. Rev. D 70, 022003 (2004) |arXiv:gr-qc/0401108| . 

M. Maggiore, Phys. Rep. 331, 283 (2000) 

P. L. Bender & D. Hils, Class. Quant. Grav. 14, 1439 (1997) 

A. J. Farmer and E. S. Phinney, Mon. Not. Roy. Astron. Soc. 346, 1197 (2003) 
|arXiv:astro-ph/0304393"] . 

D. Hils, P. L. Bender, & R. F. Webbink, ApJ 360, 75 (1990) 
R. F. Webbink, Astrophys. J. 277, 355 (1984) 

R. F. Webbink & Z. Han , in "Laser Interferometer Space Antenna: Proc. 2nd Int. LISA 
Symp.", p. 61 (1998) 



20 



[63] D. I. Kosenko & K. A. Postnov, Astron. Astrophys. 336, 786 (1998) 
[64] D. Hils & P. L. Bender, ApJ 53, 334 (2000) 

[65] G. Nelemans, L. R. Yungelson, S. F. Portegies Zwart, & F. Verbunt, Astron. Astrophys. 365, 
491 (2001) 

[66] R. Schneider, V. Ferrari, S. Matarrese & S. F. Portegies Zwart, Mon. Not. R. Astron. Soc., 



submitted, astro-ph/0002055 (2000) 
[67] G. Nelemans, F. Verbunt, L. R. Yungelson & S. F. Portegies Zwart, Astron. Astrophys., 360, 
1011 (2000) 

[68] V. Kalogera, in Gravitational Waves: Third Edoardo Amaldi Conference, ed. S. Meshkov (AIP 

Conf. Proc. Vol. CP523), 41 (2000) 
[69] D. Hils, Astrophys. J. 381, 484 (1991) 
[70] B. F. Schutz, Class. Quant. Grav. 16, A131 (1999) 



21 



VI. APPENDIX: SENSITIVITY LIMITS AND BINARY BACKGROUNDS 



This appendix summarizes for convenience the backgrounds shown in Fig. 1. More detail 
on the arguments for these numbers, and on the maximally sensitive Sagnac techniques for 
measuring the stochastic background with LISA-type interferometers, can be found in (3^. 

The approximate threshold sensitivity of the planned LISA antenna with 5 x 10 6 km arm 
lengths and for a signal-to-noise ratio S/N = 1 is shown in Figure 1. The sensitivity using 
the standard Michelson observable can be approximated by a set of power law segments: 



h rms = 1.0 x 10- 20 [f/WmHz]/Vtte, lOmHz < f (20) 
1.0 x 10- 20 /VB^, 2.8mHz < f < lOmHz 
7.8 x 10- w [(0.1mHz/ f) 2 ]/VHz, O.lmHz < f < 2.8mHz 
7.8 x 10- 18 [(0.1mi^//) 2 - 5 ]/v / Hz", OMmHz < f < O.lmHz 

where the sensitivity has been averaged over the source directions. Below 100 /iHz there is 
no adopted mission sensitivity requirement, but the listed sensitivity might be achieved for 
frequencies down to (say) 10 fiHz. More accurate modeling can be found in 

Experiments characterize backgrounds by h 2 ms , the spectral density of the gravitational 
wave strain (also sometimes denoted Sh). For cosmology, we are interested in sensitivity in 
terms of the broadband energy density of an isotropic, unpolarized, stationary background, 
whose cosmological importance is characterized by 



2 



^Cw(f) = P- 1 ^ = ^/KmsU) (21) 

where we adopt units of the critical density p c . The broadband energy density per e-folding 
of frequency, ^aw(f): is thus related to the rms strain spectral density by [57] 

hlVt GW ^ ( Kmsif) \ 2 ( f V / 22 n 
10" 8 ~ V2.82 x lO" 18 Hz-VaJ V lmHz / ' 

where ho conventionally denotes Hubble's constant in units of 100km s -1 Mpc -1 . 

Figure 1 shows the effective confusion noise from both galactic and extragalactic bina- 
ries remaining after the resolved binaries have been fitted out of the data record (see e.g. 



59l|). Essentially none of the extragalactic stellar-mass binaries can be resolved with 
LISA's sensitivity |59(] (in contrast to intense signals from an expected small number involv- 
ing massive black hole binaries). There is about a factor of three uncertainty in the estimated 
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total galactic signal level (which translates to about an order of magnitude in Q units), and 
more in the estimated extragalactic signal level. 

The normalization for various populations of galactic binaries 

(UB+WUMa+GCWDB+CV) is taken from levels estimated in |60j for the total bi- 
nary backgrounds, with estimates from 5s| for the reduction of confusion noise at higher 
frequencies by fitting out Galactic binaries. For close white dwarf binaries (C WDBs) 



60, 



61 



62, 



63, 



64, 




this value lies within a factor 2 of [62j. The GCWDBs[5<: 
include He- He, He-CO and CO-CO white dwarf binaries, as well as a few binaries containing 
the rarer O/Ne/Mg white dwarfs. 

The extragalactic XGCWDB background goes away above about 0.1 Hz, provided that 
merger-phase and ringdown radiation from coalescences are not significant. At higher fre- 
quencies, the binary background is expected to be almost entirely due to extragalactic 
neutron star binaries and 5 to 10 solar mass black hole binaries. The neutron star binary 
coalescence rate in our galaxy is taken as 1 x 10 _5 yr~ 1 . This estimate has a high uncer- 



tainty |6J3| , but gives a plausible estimate of the total gravitational wave background level, 
allowing for some additional contribution from black hole binaries. A small correction is 
included for eccentricity of the NS-NS binaries j^S], and we adopt the ratio of 0.3 between 



the extra gala ctic and galactic amplitudes 



63J. The BH-BH binaries may be the dominant 



source 1 6^ 



70]. 
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